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Abstract— The frequency domain hybrid Method of 
Moment (MoM) and Physical Optics (PO) usually requires 
a iteration between MoM and PO region. This is a 
reasonable choice for frequency domain because most of 
the computation cost is in factorizing the MoM matrix and 
while the cost of back substitution to the matrix equation 
is negligible. But it is not the case for the time domain 
version. In time domain method, the back substitution is 
the major bottle neck for computation time because it is 
done for each time step for marching-on-in-time or order 
step for marching-on-in-degree. So the iteration between 
MoM and PO region become costly because it requires to 
compute all the time or order steps for each iteration. In 
this paper, I presented a method that can avoid this 
iteration and thus speed up the computation speed for 
several times. 


Index Terms— Method of Moments (MoM), time 
domain electric field integral equation (TD-EFIE), 
physical optics (PO). 


I. INTRODUCTION 


The Method of Moment (MoM) is a widely used technique 
to solve electromagnetic field problems.[1][2] But as the size 
of the problem increases, the matrix size and the computation 
time will scale as the square and cube of the unknown 
numbers, respectively. Therefore, solution of electrically large 
problems using MoM becomes quite challenging. However, 
in many practical cases, these large objects are metallic and 
the surfaces are relatively smooth. So it is reasonable to 
hybridize a MoM methodology with asymptotic 
techniques.[3]-[8] 

In this paper, a faster time domain MoM and Physical 
Optics (PO) hybrid method is proposed. This new method 
based one the traditional theory but can solve the problem in 
a much faster manner by reducing the iterations. 


For traditional hybrid MoM-PO method, the problem of 
interest can be divided into two regions. One is the MoM 
region and the other is PO region. Assuming the PO current to 
be zero, the currents in the MoM region are solved using the 
time domain electric field integral equations (TD-EFIE), 
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which was already studied by extensive researchers[9]-[17]. 
Then the scattered field from MoM region induces a current 
on the PO region and which can be calculated by PO method. 
Hence, the PO current is known. Treating this new PO current 
as an additional source for MoM problem, the original MoM 
current can be updated. By repeating this iteration, one 
converges towards the final solution. 


But this iteration takes a significant amount of time and I 
discovered a method that can avoid this iteration but can still 
obtain the same result at no additional cost. 


In section 2, I will briefly describe the traditional time 
domain MoM-PO method. In Section 3, the proposed new 
method will be described. The cost of the new method is 
analyzed in Section 4. And finally numerical examples will be 
provided in Section 5. 


II. TRADITIONAL SOLUTION TO THE HYBRID METHOD OF 
TIME DOMAIN MOM AND PO 


For a hybrid MoM and PO problem, the structure can be 
divided into two regions. One is the MoM region governed by 
MoM formulas and the other is the PO region. For simplicity, 
we discuss only perfect electric conductor (PEC) for in this 
paper. 

For the MoM region, the boundary condition on PEC 
surface is the tangential electric field is zero. 


aaa (r, t) + E? (r, t) + poirion (r, D| = 0 (1) 


where E” (r,t), and E” (r,t), E*“’" (r,t) denote the 
scattered wave from the MoM region, the scattered wave from 
the PO region, and the excitation wave, respectively. And the 
subscript tan implies the tangential components of the fields. 
The scattered fields E(r,t) can be obtained from the current. 
Suppose we know the current in the PO region, then the 
scattered wave from the PO region E” (r,t) can be treated as 
a special form of excitation wave and the current in the MoM 


region can be solved similar to previous MoM methods. 
[9][10] 


any form is permitted without written permission by the author.* 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


If the current in the MoM region is already known, the 
current over the PEC PO region is given by the boundary 
condition of Eq. (2) and (3). In the PO region, the mutual 
interaction with PO current at other positions is neglected. The 
total magnetic field can be substituted by two times of the 
wave incident to the PO surface at the lit region Sj; and 0 in 
the shadowed region Sshadow- 


JP, 1) =x Ar, 1) (2) 
A exciatation MoM 
yr, Az 2ñx(H (r,t) +H (r, 1)) re Sj, 
0 re S shadow 
(3) 
where Heraion is the excitation magnetic field and the 


H” is the scattered magnetic field from the MoM current. 


The hybrid MoM-PO methodology can now be outlined 
using an iterative method. 


1. We first let the initial guess for the PO unknown 


current to be zero, i.e., Foy =0 


2. The MoM unknown coefficients at the k-th iteration 


J nas can be evaluated by the PO coefficients of the 


previous iteration J (it) ‘ 


3. The PO unknown coefficients at the k-th iteration J is 


can be evaluated by the MoM coefficients of the 


previous iteration J ae 


4. Ifthe results have not demonstrated convergence, then 
go back to step 2. 

In this method, we have an intention in solving the coupling 
of the coupling between MoM and PO. But in each time- 
domain MoM solution, it also involves many time step or 
order step iterations to obtain the MoM solution from an 
excitation wave. So the total iteration number will be the 
numbers of these two types of iteration multiply up. 


IHI. A METHOD TO REDUCE THE ITERATION NUMBER 


The previous method simulates the wave bounce back and 
forth in physics point of view and easy to formula. But it needs 
to solve the MoM problem for multiple times. This will not be 
a problem if this is frequency domain problem because the 
MoM matrix is the same for all iterations and thus the matrix 
factorization only need to be performed once and back 
substitution is fast. 


However, for time domain problem, the computation cost 
is different. In many cases, the major cost in solving MoM 
problem is not the matrix factorization but the backward 
substitution because it needs to be performed for each time 
step for marching-on-in-time solver or order step for 
marching-on- 
in-degree solver. So reducing the iteration number will 


significantly speed up the computation. 


If we look at the problem in a matrix format, the iteration 
between the MoM and PO coupling can be removed. As 
described in previous section, the effect of PO current is 
considered as a kind of excitation wave. 


[ Mom IL JY™ ] = [yee -MoM ] + [ 7 PO-MoM IL J” ] (4) 


in which Z“™ is MoM impedance matrix, J” the unknown 
vector of the MoM currents, J”? is the vector of the PO current 
and Z?0-MeM is the impedance matrix of the effect of PO 
currents, and Ve*<i/ation-MoM is the excitation wave. 


The PO current is calculated by the scattered wave from the 
MoM current and the excitation wave. The difference here is 
that it is not necessary to compute a matrix equation in order 
to obtain the PO current when the excitation is known. 


[ J PO | = [Vv exciatation—PO | + [ Mom -PO | [ ypitem | ( 5) 


in which VeHation-PO is the excitation wave on the MoM region 
and Z”@0™-P0 is the impedance matrix of the effect of MoM 
currents. 


In the traditional method, researchers solve Eq. (4) and (5) 
in an iterative way described in previous section. But when we 
write it in a matrix format, it is clear that Eq. (5) can be 
substitute into Eq. (4) and one can obtain the following 
equation. 


(ea 


6 
= Ba ] sf aa ieee ( ) 


The JY in Eq. (6) can be directly solved without any 
iteration. And it implies a faster simulation. After JY is 
known, J’? can be easily computed by Eq. (5). 


IV. MATRIX FILLING TECHNIQUE NOT TO INCREASE MEMORY 
USAGE 


A major difficulty in this equation in the memory space 
needed. In the traditional method, the matrix requires a storage 
space is ZY% and the size is N“°“-by-N”™, in which NY% is 
the number of unknowns in the MoM region. 


But in this formula, we need two more matrixes Z?0’™ and 
TMM EO which are N??-by-WM™ and N¥eM_by-NPO, 
respectively, in which N° is the number of unknowns in the 
PO region. The size of N”? is usually much larger than N”, 
If we need to store these two matrixes, it may require 
unaffordable amount of memory space. 


Luckily, there is a filling method that can avoid this 
expensive memory usage. We can calculate the i-th column of 
ZPO-MoM and the i-th row of Z”°”-P°, One can multiplies this 
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column to this row and add the result to the matrix ZY”, After 
doing this for all rows and columns, the multiplication of two 
matrixes of Z?OeM and ZM@eM-PO is done and no additional 
memory space is required. What is more, there are numerical 
packages such as MKL that can perform this operation in a 
very fast way. 


V. NUMERICAL EXAMPLES 


Sample numerical result is presented to illustrate the 
versatility of this method. In this work, all the structures are 
excited using a T-pulse.[18]-[22] A T-pulse is a strictly time- 
limited pulse with most of its energy concentrated in a narrow 
band. Hence the waveform by definition is causal. 


Example 1 deals with a z-polarized bicone antenna with 
reflector as shown in Figure 1. The bicone has a radius of 0.2m 
and its length is 2.1m. It also has a wire of 0.1m connecting 
the two cones in the middle. The wire has a length of 0.1m and 
radius of Imm. The excitation locates at the center of the wire. 
A 2m-by-2m square reflector is placed 2m away from the 
bicone. The MoM region contains only this bicone with 750 
unknowns and the PO region contains the reflector with 2094 
unknowns. A T-pulse with duration of 23.33 nanosecond, 
which implies a bandwidth of 300 MHz is radiated by the 
bicone. 


The current at the center of the wire is calculated with both 
traditional and proposed MoM-PO hybrid methods and 
plotted in Figure 3. It can be seen that the proposed method 
matches the traditional method well. Also the memory and 
time usage are compared in Table 1. It can be seen that the 
memory usage is not changed but the computation time is 
speeded up for 3.11 times. That is because the tradition 
method needs three iterations between MoM and PO current 
to converge for this case. 





Pe 


Figure 1. A bicone antenna with square reflector. 
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Figure 2. The excitation pulse.. 
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Figure 3. The current at the center of the wire. 
Table 1. Comparison in memory usage and computation time. 
Method Time (sec) Memory (MB) 
Traditional MoM-PO hybrid 1125 5 
method 
Proposed MoM-PO hybrid 361 5 
method 





VI. CONCLUSIONS 


A new method to compute the time domain hybrid MoM- 
PO problem is presented. This hybrid method is 
computationally efficient as it doesn’t require iterations 
between MoM and PO regions for the numerical convergence 
of the results. In all these results, a T-pulse which is strictly 
time-limited with most of its energy concentrated in a narrow 
band is used for excitation. Sample numerical result was 
presented to illustrate the versatility of this method. 
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